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Recently, many fascinating properties predicted for metamaterials (negative refraction, superlensing, 
electromagnetic cloaking,...) were experimentally demonstrated. Unfortunately, the best achievements 
have no direct translation to the optical domain, without being burdened by technological and conceptual 
difficulties. Of particular importance within the realm of optical negative-index metamaterials (NIM), is the 
issue of simultaneously achieving strong electric and magnetic responses and low associated losses. Here, 
hybrid metal- semiconductor nanowires are proposed as building blocks of optical NIMs. The metamaterial 
thus obtained, highly isotropic in the plane normal to the nanowires, presents a negative index of refraction 
in the near-infrared, with values of the real part well below — 1, and extremely low losses (an order of 
magnitude better than present optical NIMs). Tunability of the system allows to select the operating range in 
the whole telecom spectrum. The design is proven in configurations such as prisms and slabs, directly 
observing negative refraction. 

The so called metamaterials are artificial materials in which the effective medium properties, usually exotic 
and not naturally attainable, depend on the geometry of their basic constituents, rather than on their 
chemical composition 1,2 . Ranging from sub -diffraction resolution 3,4 or spontaneous emission 5 to extreme 
control over the flow of light 6,7 , many exciting and unexpected phenomena have been achieved or predicted for 
this new kind of materials. Although originally developed in electromagnetics, many of the ideas developed in this 
field have been successively extended or adapted to other ondulatory phenomena, such as acoustics, making it one 
of the most active in the engineering and physical sciences in the past few years. 

Still, however, there are many open challenges in the field. Among them, the realization of a bulk isotropic 
negative index metamaterial (NIM), with negative refraction and low losses in the optical domain 8,9 . In general 
terms, the major issue when trying to achieve such a goal is obtaining a strong diamagnetic response of the 
constituents, enough to lead to an effective negative permeability. Probably as a consequence of the success of the 
original designs operating in the microwave regime, many efforts were made to adapt them to increasingly higher 
frequencies, mainly by miniaturization 8,9 . Apart from some inherited from the original designs, such as aniso- 
tropy, many drawbacks were found in doing so, mainly related to the different behaviour of metals at optical 
frequencies, such as saturation of the magnetic response 10 or high losses associated to ohmic currents. As a 
consequence, completely different strategies were studied to obtain artificial magnetism 11,12 , as those based in 
displacement currents appearing in nanoparticle clusters due to coupling between structures 13 . However, some of 
the most successful were those which attempted to obtain it from natural magnetic resonant modes in high 
permittivity dielectrics, leading to low-loss magnetic materials 14-18 . Secondary structures or particular arrange- 
ments were needed to provide the additional electrical response, necessary to obtain doubly negative index of 
refraction 19 21 . Thus far, the maximum theoretical figure of merit (f.o.m = —Re(n e ff)/Im(n e ff)) reported is of the 
order of f.o.m.— 25 22 " 24 for metamaterials based on canonical fishnet designs; their dimensions, however, are in 
the very limit of validity of the effective medium description, and isotropy was not proven (see 25 for a review on 
this topic). 

In this work, we propose a structure that, combining electric and magnetic responses, can be used as the basic 
building block for extremely low-loss (f.o.m.— 200) isotropic two-dimensional metamaterials in the near-infra- 
red, with simultaneously negative permittivity (e) and permeability (fi) at optical frequencies. Such structure is a 
core-shell nanowire (NW) of circular cross section. Noble metals such as silver or gold can be employed to build 
the core, and high permittivity semiconductors, such as silicon or germanium, to build the shell. Artificial 
magnetism of the effective medium is achieved by exciting the lowest order Mie-like magnetic resonance in 
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the wires. We show that it is possible to tune the geometrical para- 
meters of the system to make this resonance spectrally overlap with 
the lowest electric, dipole-like, one, in turn a localized surface plas- 
mon resonance (LSPR) excitation in the metallic core. Polarization of 
the incident wave with the magnetic field along the axis of the wires 
(TE-polarization) will generate the desired responses. 

Results 

Electric and magnetic resonances of single hybrid nanowires. The 

optical properties of coated cylinders under plane wave illumination 
can, indeed, be analytically studied 26 . Under TE polarized light, the 
scattering and extinction efficiencies of these particles can be written 
as sums over different coefficients: 



2 Ekl 



-Re 



(1) 



(2) 



which can be then assigned a certain multipolar character, either 
magnetic or electric. Although the expressions are similar to those 
of a solid cylinder (i.e., without coating), the aj coefficients contain 
information about both the core and the shell (see Suplementary 
Information, Section 1). With this given polarization, the first 
zero-order coefficient, a 0 , has magnetic character, while the second 
coefficient, a u can be identified as the electric dipolar contribution in 
the 2D plane perpendicular to the cylinder axis. 

We can realize that this is the case if we study the optical properties 
of a pure silicon (Si) cylinder. Incidentally, Silicon spheres 18 have 
attracted the interest of the community, because of their potential 
use in metamaterials and interesting radiative properties. In Fig. 1 we 
have plotted the scattering efficiency of a solid Si cylinder of radius R 
= 170 nm under normal incidence of a TE-polarized plane wave. It 
can be immediately realized that the optical properties are quite 
similar to those of a pure Si sphere 18 , although the correspondence 
between Mie coefficients and electric and magnetic multipolar terms 
is not straightforward for cylinders. Thus, we need to identify the 
character of the a 0 and a x coefficients. To this end, we plot the total 
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Figure 1 | Optical properties of an infinite silicon nanocylinder. 

Scattering efficiency (black curve) of a silicon cylinder with radius 

R = 170 nm under TE-polarized plane wave incidence. Contributions of 

the two first terms in the expansion of the field are also shown. Depicted in 

the insets are the maps of the near electric field (left), and out-of-plane 

component of the magnetic field (right) at the magnetic (green, 

Iq = 1605 nm) and electric (red, X 0 = 1050 nm) resonances. The 

TE-polarized plane wave impinges from the right hand side. 



electric field and the out-of-plane component (only non-zero) of the 
magnetic field in the insets of Fig. 1 for the first two resonances. 
Similar to the case of a sphere, the magnetic resonance is the one 
excited at lower frequencies (green curve), thus corresponding to the 
a 0 coefficient. The strong circulation of the electric field is a clear 
signature of the magnetic character of the resonance. On the other 
hand, the subsequent higher-energy resonance reveals the dipolar 
electric character of the a x resonance. 

Combining a dielectric shell, which exhibits magnetic properties 
similar to those of solid cylinders, with a plasmonic core, will provide 
the additional electric resonance necessary to obtain a doubly res- 
onant structure suitable as a NIM building block (meta-atom). In 
order to make both resonances spectrally overlap, tuning of the 
geometrical parameters of the structure is needed. This approach 
has been successfully employed in the case of a core- shell sphere 27 , 
but it has to be proven for a core- shell cylinder due to the different 
behaviour of its magnetic resonance. In Fig. 2a, we have plotted the 
magnetic and electric contributions to the scattering, together with 
the total efficiency, for core-shell NWs of outer radius R out = 170 nm 
as a function of the wavelength of incidence and inner radius (Ri nt ). 

The materials used are silver (Ag) for the core, and Si for the shell. 
The electric and magnetic responses spectrally overlap over a rela- 
tively wide range, 1 0 e (1100, 1500) nm, for R int values between 60 
and 100 nm. This makes the core-shell nanowire design particularly 
robust against possible fabrication defects. Notice that, in the limit of 
^int —> 0> we recover those values of the pure Si cylinder. Although 
essentially similar, both resonances redshift in the case of Ge coating 
for the same fixed R out value, due to the slightly higher index of 
refraction of this material. The study of this case can be found in 
Section 2 of the Supplementary Information (see Supplementary 
Figure 1). 

Let us choose a particular value of R int , namely R int = 80 nm. In 
Fig. 2b the total scattering efficiency is plotted, together with the 
magnetic (a 0 ) and electric (a x ) contributions. These two contribu- 
tions are totally dominant in the frequency range in which both 
coincide. Moreover, plotting the near- field distributions of both 
the electric and magnetic fields at the frequency of coincidence of 
the resonances allows one to identify the electric and magnetic char- 
acter of the combined resonance. Indeed, a strong rotation of the 
electric field is observed, together with the signature in the magnetic 
field of a LSPR excitation in the core of the system (see insets in 
Fig. 2b). 

Metallo- dielectric core-shell NWs are, therefore, doubly resonant 
structures, having a relatively small electrical size. In the Ag@Si case 
studied, this electrical size reduces to -R 0 ut/^res ~ 0-13. These prop- 
erties allows one to envisage a possible application as building blocks 
for NIM at optical frequencies, as their spherical geometry counter- 
parts 27 . Indeed, for compact arrangements, with lattice periods of 
about the size of the diameter of the wires, the metamaterials thus 
obtained may be well described as an effective medium (depending 
also on the absolute values of the effective index of refraction, « eff , 
achieved). These compact arrangements can actually be employed in 
realistic designs, if we just impose the structures to be unconnected. 
This is a necessary requirement in order to keep the magnetic res- 
onance condition of the system unchanged. Coupling between 
neighbouring particles is small, since most of the electromagnetic 
energy is localized well inside the shell (see right inset in figure 2). 
Thus, as will be shown later, no significant coupling occurs in gen- 
eral, even if we bring the structures very close 27 . 

Physical origin of the negative index of refraction. The negative 
permeability. Before going into the study of the optical properties of 
arrangements of these core- shell nanowires, let us show that the 
doubly resonant character of these core-shell nanostructures can 
lead, indeed, to an effective negative index of refraction. Contrary 
to some other metamaterial designs in which the electromagnetic 
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Figure 2 | Optical properties of an infinite silver- silicon core-shell nanocylinder. (a) Magnetic (left) and dipolar electric (centre) contributions to the 
total scattering efficiency (right) of a 170 nm Ag covered with Si (Ag@Si) core-shell nanocylinder as a function of the inner radius, R inty and the incident 
wavelength, (b) Scattering efficiency (black curve) for a Ag@Si core-shell nanowire with external radius R out =170 nm and internal radius R int = 80 nm. 
Electric (red curve) and magnetic (green curve) contributions are also shown. Depicted in the inset are the maps of the near electric field (\E\, left), 
out-of-plane component of the magnetic field (\H Z \, center) and total electromagnetic energy density (u T , right) at the combined resonance, i.e., 
at 2 0 = 1298 nm. The TE-polarized plane wave impinges from the right hand side. 



response is such that a wave incident from free space experiences 
negative refraction in the metamaterial although the phase front 
always undergoes positive refraction 12 , the design proposed here 
is such that both the electric and magnetic material parameters are 
simultaneously negative, leading, as will be shown below, to negative 
refraction with actual phase reversal. 

Being a negative magnetization the key property to access negative 
permeablity once we neglect possible anysotropic effects, it becomes 
necessary that the basic constituents present a strong negative mag- 
netic moment: 



r x ](r)dr 



(3) 



where V represents the volume of the particle. 

If we assume the whole current to be a polarization current, 
J~(e— 1)E, the magnetic moment of each core-shell nanowire is 
simply given by: 

m=^J^ 0 [e(r)-l](rxE(r))<fr (4) 

If we now take into account the translational symmetry of the system 
along z, with the polarization current localized in the two-dimen- 
sional plane perpendicular to the nanowire's axis, we can define a 
magnetic moment per unit length, which is necessarily directed along 
the nanowire's axis m = m z z, with: 



"O 

E 

TO 



C 
QJ 
E 

o 
E 



cu 
c 

C7> 




700 



1000 



1300 
X (nm) 



1600 1900 



Figure 3 | Magnetic moment of an infinite silver- silicon core-shell 
nanocylinder. (a) The amplitude (black curves) and phase (blue curves), 
and (b) real (green curves) and imaginary (red curves) parts of the 
magnetic moment generated on a Ag@Si core-shell nanowire with 
R in = 80 nm and R out =170 nm under a TE-polarized plane wave 
incidence are depicted as solid lines. Corresponding dashed lines are the 
values for a Ag nanocylinder of radius 70 nm embedded in a Si matrix. 
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In figure 3a, we have plotted both the amplitude and phase of the 
magnetic moment generated on a Ag@Si core-shell nanowire with 
R in = 80 nm and R out = 170 nm under a TE-polarized plane wave 
incidence. Its resonant behaviour is manifested as an increase of the 
amplitude and an abrupt change in its phase. This change in the 
phase makes the real part of the magnetic moment become negative 
at frequencies above that of the resonance (1 0 = 1298 nm), see 
figure 3b. It is precisely for this frequency for which the phase 
acquires the value nil. 

The generation of such strong magnetic moment can also be seen 
as a consequence of the strong circulation of the electric displace- 
ment field in the shell of the system (see inset panel in figure 2b). To 
prove that this interpretation is correct, the magnetic moment of a 
pure Ag cylinder of radius 70 nm embedded in a Si matrix has been 
also computed (dashed lines in figure 3). This system presents a 
broad electric dipole resonance in its scattering cross section for k Q 
— 1.4 fim. Nevertheless, the magnetic properties of the sytem com- 
pletely disappear, as evidenced through its magnetic moment. 

It is easy to envisage that, when enough number of these core- shell 
nanoparticles are present in a medium, the total magnetization, 
defined as the magnetic moment per unit volume, can reach negative 
values, and so can the effective magnetic permeability of the system. 

Optical response of a periodic arrangement of doubly-resonant 
hybrid nanowires. Let us now show that it is, actually, possible to 
obtain a NIM with these nanoparticles. To do so, we arrange the 
meta-atoms in one of the most common periodic arrangements, 
namely, a hexagonal lattice. We simulate infinite slabs of different 
thicknesses to retrieve the effective index of refraction from complex 
reflection and transmission coefficients 27 " 29 . Real parts of the effec- 
tive index of refraction are, moreover, obtained by simply applying 
Snell's law. This is done under several angles of incidence (0,-). The 
retrieved parameters are then tested to fulfil Snell's law also in prism 
configurations. In the following, it will be demonstrated that the 
metamaterials obtained using Ag@Si core-shell nanowires exhibit 
not only a negative index of refraction in the optical frequency 
domain, isotropic in the plane perpendicular to the nanowire axis, 
but, moreover, extremely low losses, f.o.m.— 200 are achieved; about 
one order of magnitude better than the best reported configuration. 

Since we consider a periodic arrangement of Ag@Si core- shell 
nanocylinders in a hexagonal lattice, a highly isotropic effective res- 
ponse is expected, due to the high symmetry of the unit-cell. We take 
the geometry of the meta-atoms to be the same as in figure 2b. For 
these nanostructures there is a spectral overlap between the electric 
and magnetic resonances for k Q e (1150, 1500) nm. We expect a 
metamaterial with these basic building blocks to present a left- 
handed (negative-index) frequency window in which propagation 
is allowed. To obtain this effective response, relatively high filling 
fractions (/) are expected to be needed 27 . Therefore, basic constitu- 
ents are brought close together, fixing the surface-to-surface distance 
of every two neighbouring nanoparticles to be d s _ s =10 nm. For a 
hexagonal lattice this implies a filling fraction off — 0.856. We have 
performed full numerical simulations using a Finite Element Method 
(FEM) commercial software (COMSOL Multiphysics 4.2). Slabs are 
considered infinite along transverse direction by imposing appropri- 
ate periodic boundary conditions to the system (see section Methods 
for further details). Let us now define a basic layer of these slabs to be 
such that different thicknesses are obtained by adding a different 
number of these layers, without need of any lateral translations. 
They are infinite in the transverse directions, and we call them unit 
constituents of the slabs (u.c). 

In Fig. 4a, the obtained transmission spectra under normal in- 
cidence illumination are depicted. Results for slabs of several 



thicknesses, expressed in terms of the number of unit constituents, 
N, are included, the basic u.c. being described in the inset of figure 4b. 
A frequency window in which propagation is allowed is present, 
which coincides with the range in which both resonances overlap. 
Thus, coupling between neighbouring nanostructures does not seem 
to have a high impact, since the resonances are not critically modi- 
fied. Recall that, although d s _ s is small compared to the size of the 
structures, the highest near field intensities are located in the surface 
of the core. Since core-to-core distances, d c _ c , are of the order of d c _ c 
~ 2 (i^ — R int ), coupling between plasmon modes is very weak, and 
the resonances are not disturbed. Modulations of the intensity inside 
the transmission window correspond to well known Fabry-Perot 
modes inside the slab, due to multiple reflections. The transmission 
peak arising at higher energies (2 ~ 1050 nm), slightly redshifted 
with respect to the electric-quadrupole-like resonance, almost com- 
pletely disappears for sufficiently thick slabs, as evidenced in the 60 
u.c.-thick slab case (see the inset in Fig. 4a). The obtained values 
within the transmission frequency range decrease with the thickness 
of the slab due to absorption in the system. Nevertheless, there is still 
a significant amount of energy transmitted through, even for such a 
thick slab as the one with 60 u.c. (i.e. for propagation lengths inside 
the NIM of about 36 fim). Within this frequency window, negative 
refraction is thus expected with low associated losses. Before retriev- 
ing the values of the effective parameters, let us show that this NIM is, 
moreover, highly isotropic. To do so, we have performed the same 
study for several angles of illumination, 0 t . Results obtained for the 60 
u.c.-thick slab and 0 t e (0°, 45°) are shown in Fig. 4b. For low angles of 
incidence, 0,- e (0°, 35°), transmission curves are almost identical. For 
increasingly higher angles of incidence, Q { e (35°, 45°) a Fabry-Perot 
modulation of the transmission is observed in the high frequency 
limit. 

Negative refractive index retrieval. As previously mentioned, we 
took three alternative approaches to determine, within the trans- 
mission window, Re(n e ff) from simulations. The results thus ob- 
tained are depicted in the lower panel of Fig. 5. First, we have 
applied standard retrieval procedures 28 ' 29 , which allows one to deter- 
mine the effective parameters from complex reflection and trans- 
mission coefficients (or S-parameters) under normal incidence 
illumination, assuming an equivalent homogeneous slab. The blue 
dashed curve represents the obtained values. We have then applied 
Snell's law to those slab configurations in which the angle of 
incidence is 6j 0. The results for the particular case of a 20 u.c.- 
thick slab under 0 t = 20°, 30° illumination are plotted as red and 
green squares. Additionally, we simulated a rectangular prism with 
an angle 0 pr i Sm = 30° and applied Snell's law: the thus retrieved values 
are indicated by black circles. Both insets represent the electric field 
component parallel to the first interface in the slab and prism 
configurations for X ~ 1215 nm. Negative refraction with extre- 
mely low losses is clearly evident. Moreover, phase reversal, charac- 
teristic of systems having simultaneously e < 0 and fi < 0, can be 
observed by plotting dynamically this component in a full harmonic 
cycle (see Supplementary Videos 1 and 2). Errors in those values, 
retrieved graphically from Snell's law, have been computed by 
assuming a systematic error in the determination of the angle of 
±2°. It is apparent from Fig. 5 that the retrieved values are 
consistent regardless of the angle of incidence, indicating, as 
expected, that the obtained metamaterial is highly isotropic. 

It is important to mention that, although the standard retrieval 
procedure based on S-parameters allows one to compute the index of 
refraction in the whole frequency range, the effective medium 
description of the system breaks down at longer wavelengths inside 
the transmission window. If somehow counter- intuitive, this effect is 
a consequence of the fact that dRe(n e ff(co))/dco > 0 (non-anomalous 
dispersion), leading to higher absolute values at longer wavelengths, 
which result in shorter effective wavelengths inside the metamaterial, 
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Figure 4 | Transmission curves through metamaterial slabs made of a regular hexagonal arrangement of Ag@Si core-shell nanowire with external 
radius .R out = 170 nm and internal radius R int = 80 nm. Centre-to-centre distance between neighbours is L = 350 nm. (a) Transmission curves for 
several thicknesses under normal incidence illumination, corresponding to different symbols. The inset shows the case of a slab with a thickness of 60 unit 
constituents (u.c). (b) Transmission curves through a 60 u.c- thick slab under several angles of incidence. 
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Figure 5 | Effective index of refraction of a metamaterial built with Ag@ 
Si core-shell nanowires with .R int = 80 nm and .R out =170 nm arranged 
in a hexagonal lattice. Lower Panel: Retrieved real part of the index of 
refraction. The blue discontinuous curve represents the values obtained 
from complex reflection and transmission coefficients under normal 
incidence. Black circles represent the retrieved values from Snell's law in 
the prism configuration. The incidence angle is fixed, 6i = 6 prism = 30°. The 
corresponding inset shows the y-component of the electric field 
(only non-zero of the incident wave) for the particular case X ~ 1215 nm. 
Red and green squares represent the values obtained from Snell's law in the 
slab configuration. Angles of incidence are Q { = 20° (red) and Q { = 30° 
(green). The corresponding inset represents the x-component (parallel to 
the first interface) of the electric field for X ~ 1215 nm. Error bars 
computed assuming a systematic error in the determination of the 
refracted angle of ±2°. Upper Panel: Imaginary part of the index of 
refraction obtained from complex reflection and transmission coefficients 
(black curve) and the corresponding figure of merit, f.o.m., as defined in 
the text (blue curve). 



A e ff. In this way, at k 0 = 1 150 nm, the effective wavelength inside the 
structure can reach 2 eff = X 0 /\Re(n e ff)\ — 15 /nm, i.e., up to 40 times 
the characteristic period length within the structure. By contrast, at 
X 0 = 1500 nm, 2 eff = A 0 /\Re(n e ff)\ — 575 nm, so that the effective 
wavelength inside the system is smaller than two times the period. 
Therefore, when the retrieved refractive index becomes very large 
(despite being negative), the effective wavelength inside the metama- 
terial becomes correspondingly small, thus becoming of the order of 
the lattice constant. At such limit, we do not fully trust the resulting 
refractive index, at least from the point of view of proper NIM beha- 
vior, thus establishing a cutoff in Fig. 4. 

Losses and Figure of merit. With respect to losses, there is a strong 
evidence that the metamaterial is very weakly absorbing, as deduced 
from the fact that the transmission values are as high as 0. 1 within the 
negative-index band, with about the same portion of reflection, for 
propagation lengths inside the metamaterial of about 36 jdm. This 
extremely low attenuation of waves propagating inside the NIM is 
also evident in the electric near field maps (see insets in figure 5). One 
can actually compute the expected values of the imaginary part of the 
index through the standard S-parameter retrieval procedure: The 
results are plotted in the upper panel of figure 5. While values as 
extraordinarily low as 7 X 10" 3 are achieved within the transmission 
window, out of it, the metamaterial essentially behaves as a metal, 
with Re(n e ff) < Im(n e ff), reflecting almost entirely the incident wave. 
Inside the negative-index band, values of the figure of merit (f.o.m = 
—Re(n e ff)/Im(n e ff)) can, in turn, be as high as f.o.m.— 200. In the 
particular case of Re(n e ff) = — 1, f.o.m.— 85, i.e., 7m(n e ff) = 0.012. 
Recall that the f.o.m. in the case of core-shell nanospheres is 
substantially lower 27 , stemming from their much larger absorption 
efficiency, as compared to that of core-shell cylinders. 

To further test the validity of this proposal, similar analysis has 
been performed to the case of Ag@Ge core-shell nanowires arranged, 
instead, in a square lattice (keeping the d s _ s constant). Results indi- 
cate that the so obtained metamaterial behaves as a NIM in the 



SCIENTIFIC REPORTS | 3 : 1507 | DOI: 1 0.1 038/srep0 1507 



5 



1 0 e (1.3, 1.45) fim frequency range. Values for the real part of the 
index are slightly lower and associated losses higher (see Section 2 in 
the Supplementary Information). Concerning the metal core, any 
other metal with a good 'plasmonic' behavior can be used, like Au 
or Cu, leading indeed to a quite similar low-loss NIM behavior. 

Discussion 

As showed here, metallo- dielectric core-shell nanowires of circular 
cross section, are very promising candidates as building blocks of 
highly isotropic, negative index metamaterials operating at optical 
frequencies, with extremely low losses exhibiting figures of merit 
f.o.m.— 200 (an order of magnitude larger than previously proposed 
NIMs). Although here demonstrated to operate at X 0 e (1.15, 
1.35) fim for Ag@Si with R int = 80 nmand.Rout =170 nm arranged 
in a hexagonal lattice, and at 1 0 e (1.3, 1.45) fim for Ag@Ge nano- 
wires with R int = 90 nm and R out = 170 nm arranged in a square 
lattice, negative refraction in the whole telecom frequency range can 
be achieved by tuning the geometrical parameters. We have also 
verified that the low-loss NIM behavior can be scaled down to 
800 nm by introducing a thin layer of a low index material (i.e. silica) 
in between the Ag core and the Si shell. Theoretically, such thin layer 
would produce a slight blueshift of the dipolar electric resonance; 
correspondingly, the magnetic resonance could be properly matched 
by tuning the shell dimensions. To further downscale the design into 
the visible, a nearly lossless semiconductor should be employed. 
Moreover, the operating physical principles can be applied to obtain 
negative refraction in completely different ranges of the electromag- 
netic spectrum (IR and THz) by using appropriate materials. 

We also anticipate that the proposed design, though challenging, is 
no doubt amenable to fabrication at present. For instance, nanopor- 
ous silicon templates can be easily fabricated in arrays with long- 
range order in various manners 31,32 ; indeed, core-shell Si NWs 33 or 
even hollow Si nanowires 34 have been very recently obtained with 
suitable dimensions. For our proposed NIM design, filling those 
holes with Ag or Au would then suffice. Recall that many other 
high-index, nearly lossless semiconductors exist that will lead to 
the same NIM behavior. In this regard, for instance, hollow AllnP 
semiconducting nanowires have been patterned by metal- organic 
vapor phase epitaxy in a photoresist matrix 35 ; filling the hollows with 
a noble metal prior to photoresist removal would complete the pro- 
cess. On the other hand, it should be mentioned that, if needed for 
fabrication purposes, a thin shell of silica (or any other low-index 
dielectric, as mentioned above) could be deposited as an intermediate 
layer between the noble metal and the semiconductor shell. 

In addition, semiconducting core-shell nanowires have been fab- 
ricated by a two-step growth process combining selective-area and 
vapor-liquid-solid epitaxy 35 . Furthermore, we anticipate the robust- 
ness of the proposed design against NW size imhomogeneities and 
array defects/disorder; this stems from the broad spectral width of 
the NW resonances and resulting negative-index bands, along with 
the fact that the latter does not critically depend on the NW particu- 
lar arrangement. 

Methods 

The scattering properties of solid and core- shell cylinders were computed analytically 
following references 26 30 (see also Supplementary Information, Section 1). 
Corresponding near field maps were also computed analytically with this formalism. 

Full electromagnetic simulations of propagation of waves inside slabs and prisms 
were carried out using a Finite Element Method commercial software (COMSOL 
Multiphysics 4.2). In the case of slabs, Bloch boundary conditions are applied at the 
lateral ends of the simulation domain to reproduce an infinite system. At the 
boundaries ahead and behind the slab, complex transmission and reflection coeffi- 
cients (S-parameters) are computed. We use the built-in port boundary condition to 
this end. The effective material parameters are, then, retrieved from the S-parameters 
according to 27 ' 29 . Although the imaginary part of the effective index quickly converges 
with increasing number of functional layers, retrieval of the real part implies an 
ambiguity in the election of the branch of the complex logarithm 29 . To get rid of such 
ambiguity several slab thicknesses were simulated. In particular, we simulated slabs 
made of 2, 3, 4, 5, 6, 7, 10, 20, 30, 40 and 60 functional layers (not all shown in figure 4). 



Silicon was assumed lossless in the frequency range studied, with a refractive index 
of n S i = 3.5. Nonetheless, to make sure that the small imaginary part of Si is not 
affecting the resulting f.o.m., numerical calculations including it were carried out, 
yielding the same f.o.m. See Supplementary Information for a comparison of the 
optical properties of core- shell nanowires when losses are taken into account, as well 
as the retrieved imaginary part and figure of merit. Silver material properties were 
taken from experimental values from reference 36 . 
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